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ABSTRACT 

Reservoirs produced by a mechanism termed water drive. Often this is called natural water drive to 
distinguish it from artificial water drive that involves the injection of water into the formation. 
Hydrocarbon production from the reservoir and the subsequent pressure drop prompt a response from 
the aquifer to offset the pressure decline. This response comes in a form of water influx. This paper is to 
determine water influx in reservoir in Niger Delta using van everdingen-hurst unsteady-state model . In 
doing this, the theoretical aspect of these models were looked at and validation of production data was 
carried out. .then the application of the model was also considered by using production data . The 
results show that estimated value of water influx was 6305723 1.92bbl for one year and 2250770727 
bbl at the end of seven years. The results were also compared with the company values. It is clear from 
the graph there is little disparity between these models. From the t-test carried out we concluded that 
there is significant difference between the two models. For exact solution Van Everdingen-Hurst is 
suitable. However, for easy and approximate solution, Company(Hurst-VaneEverdinge-Odeh) Models 
is preferred. 
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INTRODUCTION 

There are several models used in Petroleum industry to estimate water influx, this include the 

following(Tarek.HAhmed:Reservoir Engineering Handbook): 

*Pot Aquifer 

*Schilthuis' steady- state 

*Hurst's Modified steady-state 

*The Van Everdingnen-Hurst Unsteady-state 

*The carter-Tracy Unsteady-state 

*Fetkovich's Method. 

In this paper, water influx in reservoir will be determined using van everdingen-hurst unsteady-state model. The 
values obtained will be compare with the company model This particular paper is particular seeks to empower 
the reservoir Engineers and those involve in reservoir management to know suitability of the model. 

RESEARCH METHODOLOGY 

In this section, the focus is on the design of the research, area of study and method of data collection and method 
of analysis. 

RESEARCH DESIGN 

In this work, field "X" was used as a case study and application of the model will be demonstrated .The 

following steps will be taken in using the models: 

Step 1, Determine water influx constant B 

Step2, For each time period, calculate AP 

Step3, Calculate dimensionless times t D 
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Step4, For each t D computed in step3, calculate a dimensionless cumulative water influx 

Step5, Calculate W e 

Step6, Summary of the final results 

AREA OF STUDY The work is focused on Van Everdingen-Hurst Model . 

METHOD OF DATA COLLECTION 
The data for this work was obtained from record. 

METHOD OF DATA ANALYSIS 

In analyzing the data collected, graph of reservoir pressure and time was drawn and determination of drive index 
was carried out. 

DETERMINATION OF DRIVE INDEXES 

Oil is produced from many reservoirs by more than one natural production mechanism .In managing such a 
reservoir; it is useful to make an estimate of the fraction of total production that should be attributed to each of 
the mechanisms that are present. Pirson has proposed the use of drive indexes to answer that question. The 
derivation is based on the relationship; 

(oil expansion) +(free gas expansion) + (net water influx) + (net injection) = (Production) 

If this equation is written in mathematical form,then divided by the production term on the right side of the 
equation relationship is obtained in the form (L.F.Koederitz,A.H. Harvey and M.Honarpour ,1989): 

DDI+SDI+WDI+IDI=0 1 

Where the solution gas drive index (also called the depletion drive index or DDI) may be expressed by the term 

NJBo+tesi-Bjql-NBrt 2 

The gas cap expansion drive index(also called the segregation drive index SDI) may be written: 

3 

GJ3# — GB GI 

The water drive index(or WDI): 

w e -w p B v , 4 

N p E +G p E g -N P R s B g 

The injection drive index (IDI) is: 

I-/ 1 B lt .+G;B Jj 5 

RESULTS AND ANALYSIS 

APPLICATION AND RESULTS OF VAN EVERDINGEN-HURST MODEL 

The application of this model will be illustrated using the information in Tables 1.0 and Table2.0 .Tablel.O 

gives information on estimated properties of the aquifer.Table2.0 gives information on Pressure History at the 
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Reservoir/Aquifer Boundary (Dake,L.P.,1978).As reservoir fluids are produced, reservoir pressure 
declines. Figl.O illustrate reservoir pressure declination. 

Tablel.O The estimated properties of the Aquifer. 






M(C P ) 





K(md) 


QCpsia 1 ) 


h(ft) 


r r (ft) 


0.23 


0.1728 


160 


2180 


7xl0" 6 


212 


6561.6797 



Table2.0:Pressure History at the Reservoir/ Aquifer Boundary. 



Pressure(Psia) 


3452.30 


3328.8 


3263.91 


3142.93 


2999.57 


2870.8 


2758.08 


2714.24 


Time(days) 





365 


728 


1090 


1455 


1820 


2185 


2551 


W e (MMbbl) 





1.52 


4.18 


8.13 


13.0 


18.9 


25.75 


32.52 



4000 

3UUL) 

Press(Psia) 7( ™ 

1000 









♦ 3 ressure{Ps s) 
near {Pressure{P£iai) 



moo ?ono moo 

Time{Dayl 



Figl.O Reservoir pressure decline 

DETERMINATION OF DRIVE MECHANISMS AND DRIVE INDICES 
The drive indices are calculated as follow. 

For solution gas drive index (also called the depletion drive index or DDI) may be expressed by the term: 



2.63frlQ e [1.7165 + (1378.99 - 1266.3l)8.1tolQ- 4 ] - 2,63x1Q 8 :e1.771B 
5.16^10^1.7165 +7.71.t10"^S.19.t10-* - 5.16%lCxl266.31*S.191C-* 

_9.73jtlO = QQ99]^ 
9.82 10 7 

For gas cap expansion drive index(also called the segregation drive index SDI) may be written 

GB 3 - GBcj 



.%B C ^B B -N P R S B 3 

1.1S?c1Q :l2 ?cS.19^10-' 1 ' 1.1S^10 12 ?c7.6?c10-' 1 ' 

5,ieA-irj 7 'j:1.7165 +- 7.71xl0 1D ;cS,19;cltr 4 - 5,16x-irj 7 jcl2&e.31;cB.19irj " 



6.30xlG 7 
9.82xl0 7 " 



: 0.642 
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For water drive index(or WDI): 



W 9 — W V B W 



3.25xl0 7 - 8.80a:10 5 x1.03 

5,16xia 7 jl,7165 +- 7,71;clQ 1D jr3,19;clO~ 4 - 5,16;eia 7 .cl25&, 31*8,1910" 

3.16xl0 7 

= 0.322 



9.82xl0 7 
For injection drive index (IDI) is: 
WjlAv -t- gjrJTg 

0x1.03 + 0^8. 19xl0 -4 



.=0 



Therefore from equation 3.0 

0.0991+ 0.642 + 0.322 + = 1.06 

The calculation of drive indices show that production of this reservoir involve combination of drive 
mechanisms and also show that there was water influx. 

Using the information in Tables 1.0 and 2.0. 

Stepl, Determine water influx constant B use equation 10.0 

B = 1.1190C t hr 1 2 ( i ^) 10.0 

= 1.119(0.23)(7xl0- 6 )(212)(6561.6797) 2 ^— j 

=7308.687RB/Psi 

Step2, For each time period, calculate AP using Eq.l 1.0 

^fOV.-z-JNJ 1L0 

For n= 1 

APi = tC3452.3-3328.8J = 61.75 
For n=7 
AP 7 = -C2S7Q.S - 2714,2*) = 78,28 
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Step3, Calculate dimensionless times t D using Eql2.0. 



: Q.QQ<j33fVTf _ 



-12.0 



0.00633(2180) £ 



= 1.1466t 



For example at t = 365days 



t D = (1.1466)(365) = 418.51 



Step4, For each t D computed in step3, calculate a dimensionless cumulative water influx. Because we are 
assuming an infinite-acting aquifer,we can use equations 13.0 since t D is greater than 0.01 and less than 
200(John,L. and Robert. A. W, 1996) . The value of t D determines which equation to use. For example 
at t =365days(n=l),t D =418.51 .Since t D >200, 

so we use Eq.13.0 

Qm',t D ) = — = 139./ 2 

Table3.0 summarizes the intermediate results from step2 through 4 

Table3.0 Summary of intermediate results using Van Everdingen-Hurst Model 



N 


Time(days) 


t D 


P (Psia) 


AP n (Psi) 


Qp D (tD„) 











3452.30 








1 


365 


418.51 


3328.8 


61.75 


139.72 


2 


728 


834.72 


3263.91 


94.20 


250.71 


3 


1090 


1249.79 


3142.93 


92.94 


354.43 


4 


1455 


1668.33 


2999.57 


132.17 


454.9 


5 


1820 


2086.81 


2870.8 


136.65 


552.49 


6 


2185 


2505.32 


2758.08 


120.75 


647.91 


7 


2551 


2924.98 


2697.49 


78.28 


741.85 



Step5, Calculate W e using Eq.14.0 

W e (t Dn ) = B%if= 1 APiQl>D(Vti-l)D 

For example, at n=l 

W e (t D1 ) = B£j =1 APiQpDft-ti.OD 

=B[AP 1 Q PD (t 1 -t ) D = B[AP, Q PD (t D1 )=7308.687(61.75x6.3269) 
=6305723 1.92bbl 
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Similarly, for n= 6W e (t D6 ) = B^f- x h? & m {t 6 -t iA ) D 150 = B[AP 1 Q PD (t 6 -t () ) D + AP 2 Q PD (t 5 ) D + AP 3 Q PD (t D4 ) 



+ AP 4 Q PD (t D3 ) + AP 5 Q PD (t D2 ) + AP s Q PD (t D1 )] 

=7308.687[(61.75)( 647.91)+(94.195)( 552.49)+(92.935)( 454.9)+(132.17)( 354.43)+(136.65)( 
250.71)+(120.745)( 139.72) ] 

=1 697860 109bbl 

Table4.0 summarizes final results using Van Everdingen-Hurst Model. 

Table4.0 Su mmary of final results using Van Everdingen-Hurst Model 



N 


Time(days) 


P (Psia) 


W e (bbl) 








3452.30 





1 


365 


3328.8 


63057231.92 


2 


728 


3263.91 


209342476.9 


3 


1090 


3142.93 


427474222.2 


4 


1455 


2999.57 


754586564.6 


5 


1820 


2870.8 


1185013265 


6 


2185 


2758.08 


1697860109 


7 


2551 


2697.49 


2250770727 



Table5 .0 Van Everdingen Hurst and (Everdingen&Odeh(MMBBL))Company Models 



N 


VanEverdingenHurst(BBL) 


Everdingen&Odeh(MMB B L) 











1 


63057231.92 


1.52 


2 


209342476.9 


4.18 


3 


427474222.2 


8.13 


4 


754586564.6 


13 


5 


1185013265 


18.9 


6 


1697860109 


25.75 


7 


2250770727 


32.52 



Z.5_t09 Fig 2 .0 Comparison between Van Everdingen Hurst and CormanvfHurst 



We 




Everdingen 
&0deh (BBL) 



Fig2.0 shows comparison between Van Everdingen Hurst and Company Models 
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DISCUSSION 

.In analyzing the data used for the application of this models, reservoir pressure was plotted against time as 
shown in fig 1.0. Fig 1.0, illustrates the decline in the boundary pressure as a function of time for a radial 
reservoir -aquifer system. If the boundary pressure is reduced, pressure drop will be imposed across the aquifer. 
Water will continue to expand and the new reduced pressure will continue to move outward into the aquifer. 
Since these pressure waves are assumed to occur at different times, they are entirely independent of each other. 
Thus, water expansion will continue to take place as a result of the first pressure drop, even through additional 
water influx is also taking place as a result of one or more later pressure drops. This is essentially an application 
of the principle of superposition (Craft,B., Hawkins. M and Terry. R ,1991). In order to determine the total water 
influx into a reservoir at any given time, it is necessary to determine the water influx as a result of each 
successive pressure drop that has been imposed on the reservoir and aquifer. In calculating cumulative water 
influx it is necessary to calculate the total water influx from the beginning. This is required because of the 
different times during which the various pressure drops have been effective 1 . Van Everdingen-Hurst 
methodology provides the exact solution to the radial diffusivity equation and therefore considered the correct 
technique. 

The final results for the model is shown in table4.0 above. The difference between the Company(Hurst- 
VaneEverdingen-Odeh) technique and the Van Everdingen -Hurst technique is shown in table5.0. Graphically, 
fig2.0 compare the model with the company model(Hurst-van Everdingen-Odeh) . It can be seen that there is 
little disparity (AndyJ.J., 2000) between the two models 

CONCLUSION 

In concluding this paper it is once again necessary to highlight what the original aim of the paper was in the first 
place. The aim of the work is to determine water influx using Van Everdingen-Hurst model. Objective of this is 
to determine suitability of the model in estimating water influx in Niger Delta .This will assist the reservoir 
Engineers and production managers in during reservoir management. This aim to a large extent could be said to 
have been achieved. 

RECOMMENDATION 

What has been presented in this paper is perhaps a first attempt at providing the reservoir Engineers and 
production managers a means of determining which of the model to use when calculating water influx and as 
can not be said to be conclusive. Consequently, the following recommendation for further research work have 
been made. 

1 The model can be compare with Carter-Tracy method 

2 I also recommend that the model should be compared with fetkovich model. 

NOMENCLATURE 

C, = Total aquifer compressibility. Psi" 1 

C f = Aquifer rock compressibility Psi" 1 

C w = Aquifer water compressibility Psi" 1 

h =Net formation thickness, ft 

e w = water influx rate from aquifer, bbl/day 

P aq = Aquifer pressure, psia 

P D = Dimensionless pressure 

P r = Pressure at aquifer/reservoir interface. Psia 

AP = Dimensionless difference between initial aquifer pressure and pressure at original reservoir/aquifer 

boundary, Psia. 

r r = Radius to aquifer/reservoir interface 

r a = Radius of the aquifer 

<D= Porosity of the aquifer 

W e = Cumulative water influx. RB 

Z D = Dimensionless vertical distance 

K v = vertical permeability 

K h = horizontal permeability 
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W ei = Initial encroachment volume in a aquifer, RB 
W; = Initial volume of water in a aquifer,RB 
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